
 *************************************************************************************
 *************************************************************************************
 **                                                                                 **
 **                 Welcome to Point-by-Point Regression (PPR_v2c3)                 **
 **                                                                                 **
 **   A program for producing fully nested reconstructions of gridded climate data  **
 **   with parametric prediction intervals or optionally (if you have enough time)  **
 **   semi-parametric prediction intervlas based on the maximum entropy bootstrap   **
 **                                                                                 **
 **         Programmed by:  Edward R. Cook, Lamont-Doherty Earth Observatory        **
 **        My answer to the spatial reconstruction of climate from tree rings       **
 **                  Version 2.0b (well not really):   12/23/2012                   **
 **                                                                                 **
 **              The program dimensions are currently set as follows:               **
 **                                                                                 **
 **         mxg  (maximum number of grid coordinants to read in)     =  1000        **
 **         mxp  (maximum number of climate data years to read in)   =   150        **
 **         mxy  (maximum number of years of tree-ring data to use)  =  2100        **
 **         mxs  (maximum number of tree-ring series x their lags)   =   500        **
 **        mxtst (maximum number of ensemble members to produce)     =    30        **
 **         mxb  (maximum number of MEbootstrap replications)        =   300        **
 **                                                                                 **
 *************************************************************************************
 *************************************************************************************

 Okay, have you read this? Is your recon run the right size to fit? <y>/n: 

 *************************************************************************************
 *                                                                                   *
 *  Before proceeding now, have you constructed your tree-ring headers correctly as  *
 *  the program requires? Here is what a typical tree-ring header MUST look like:    *
 *                                                                                   *
 *  PK    CHP1JU   JUEX   1640  1993   36 20    74 02   CHAPROT CHAP1 ESPER JUNIPER  *
 *                                                                                   *
 *  where:                                                                           *
 *         PK   is the location code of the chronology (e.g. country, state)         *
 *       CHP1JU is the tree-ring chronology site id code                             *
 *        JUEX  is the four letter tree species abbreviation                         *
 *        1640  is the first year of the chronology                                  *
 *        1993  is the last year of the chronology                                   *
 *       36 20  is the latitude in degrees and minutes of the site location          *
 *       74 02  is the longitude in degrees and minutes of the site location         *
 *      CHAPROT ... is the tree-ring chronology site name (max 47 characters)        *
 *                                                                                   *
 *  If the answer is NO, you must do this before running the program by creating a   *
 *  list file of all headers of all chronologies organized by geographic location    *
 *  and these headers must be identical to the headers of the chronologies read in.  *
 *  The program will check the headers in the list file for spacing and read format  *
 *  errors. If any are found, the program will not run until you fix those errors.   *
 *                                                                                   *
 *************************************************************************************

 Okay, have you carefully read this stuff above? If you are a NEWBIE, you better!

 So, create tree-ring headers now? y/<n>         ==> 

 Use an existing run time parameters file? <y>/n ==> 

 Enter run time paramters file name (/ to exit)  ==> PPR_v2c3_Params_File.txt

 Run time parameters file opened                 ==> PPR_v2c3_Params_File.txt                                    

 Climate variable file opened                    ==> Partab_Monthly_Flows.txt                
 Climate variable format to use                  ==> (i4,12f12.2)                            

 Climate variable lat-lon file opened            ==> Partab_flowXY.txt                       
 Lat-Lon read format to use                      ==> (f10.3,f9.3,a)                
 Lat-Lon coordinates order read                  ==> latlon
 Total number of climate grid points read        ==>     1
 Number of climate grid points to process        ==>     1

 Tree-ring header file opened                    ==> PK_crns_list.txt                        
 Number of tree-ring header formats checked      ==>   39
 Number of tree-ring headers formally read       ==>   39

 Calibration period first and last years         ==>  1975  2004
 Verification period first and last years        ==>  1962  1974

 JFMAMJJASONDJFMAMJJASOND                
 000011113               

 tmo,ism,iem,nmo:    9   5   9   5  5.0

 Season to be reconstructed: Months  5  9  M5S

 t=1 regression model order selection criterion chosen

 No backtracking to best verified model used.

 No adaptive power transformation applied to predictands.

              DVAR  MON   DIST  MINS  MINFY  MINC  ARMOD  NLAG  SCREEN  NTAIL  SCALE  SAVE    WGT   NBOOT
 Test:   1    FLOW  M5S    500    1   1800     5     0      2     90      2      2      2   0.000     300

   NUM   GPT    LAT    LON   MON  MOD  IFYR  ILYR   DIS  NSR  NTR  EIG  EV1  ORD   RSQ    CVRE   NYV    R2     RR      RE      CE       RMSE
     1     1  35.72   74.63  M5S    1  1800  2004   291   48   15   15    4    3  0.554   0.494   13  0.449  0.445   0.431   0.406   363.626
 
 *** Exit PPR_v2c3 ***

